Origin of coherent magnetic fields in high redshift objects 
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ABSTRACT 



Large scale strong magnetic fields in galaxies are generally thought to have 
been generated by a mean field dynamo. In order to have generated the fields 
observed, the dynamo would have had to have operated for a sufficiently long 
period of time. However, magnetic fields of similar intensities and scales to the 
one in our galaxy, are observed in high redshift galaxies, where a mean field 
dynamo would not have had time to produce the observed fields. Instead of a 
mean field dynamo, we study the emergence of strong large scale magnetic fields 
in the first objects formed in the universe due to the action of a turbulent, helical 
stochastic dynamo, for redshifts 5 < ^ < 10. Ambipolar drift plays an important 
role in this process due to the low level of ionization of the gas, allowing a large 
scale stochastic dynamo to operate. We take into account the uncertainties in the 
physics of high redshift objects by examining a range of values for the parameters 
that characterize the turbulent plasma. By numerically integrating the nonlinear 
evolution equations for magnetic field correlations, we show that for reasonable 
values of the parameters in the time interval considered, fields can grow to high 
intensities (~ 10~^ G), with large coherence lengths (~ 2 — 6 kpc), essentially 
independent of the initial values of the magnetic field. 

Subject headings: magnetic fields, high redshift objects, turbulence 
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1. INTRODUCTION 

Magnetic fields have been observed in all known structures of our universe, from the 
Earth to superclusters of galaxies, spanning a wide range of intensities from ~ /xG in 
galaxies and galaxy clusters to ~ 10 G in neutron stars. The origin of these fields in 
large structures, such as galaxies and clusters of galaxies, remains an unsolved problem. 
Many physical processes have been proposed to explain the origin and evolution of these 
fields (see reviews Grasso & Rubinstein 2001; Widrow 2002). The processes suggested can 
be divided into two main classes: 1) cosmological mechanisms; and 2) local astrophysical 
processes. Until now, none of them has provided a satisfactory explanation for the 
generation of the magnetic fields. 

A mean field dynamo is commonly invoked to explain the fields observed in our galaxy 
and in small redshift galaxies (e.g., Zel'dovich et al. 1983; Moffat 1978). In order to have 
attained the observed intensities, the dynamo would have had to have operated for a time 
on the order of the age of the universe. However, the presence of equally intense and 
coherent fields in high redshift galaxies (Carrilli & Taylor 2002), where the mean field 
dynamo would not have had enough time to amplify the field to the observed values, casts 
doubt on the mean field dynamo paradigm as the preferred generation mechanism. Fields 
of similar intensity and coherence to those in the Milky Way, have also been detected in 
high redshift damped Lyman — a. systems (Wolfe et al. 1992). 

In this paper, we are concerned with the origin of strong, large coherent magnetic fields 
in high redshift structures. We relate the origin of magnetic fields to the physical conditions 
in the early universe: 1) a dense low ionized plasma; and 2) appreciable turbulence due to 
the observed high star formation rate (Lanzetta et al. 2002). 

The formation of the first stars and quasars marked the beginning of the transformation 
of the universe from a smooth initial state to its present lumpy state. In the bottom-up 



-4- 



hierarchy of cold dark matter (CDM) cosmologies, the first gaseous clouds collapsed at 
redshifts 2; > 10 and, subsequently, fragmented into stars due to molecular hydrogen cooling 
(Barkana & Loeb 2001). These collapsing objects then fragmented into many clumps, which 
had typical masses of ~ 10^ — 10 ^M©. Very massive stars have hfetimes of ~ 3 x 10^ years 
and end their lives as supernovae. 

Recently, Lanzetta et al. (2002) showed that the incidence of the highest intensity star 
formation regions increases monotonically with redshift. Their observations indicate that 
star formation in the early universe occurred at a much higher rate than was previously 
believed. Therefore, the rate of occurrence of supernovae would have also been much higher 
in the past than at the present. Supernovae shocks disturb the plasma in which they 
are immersed, producing turbulent motions of the gas. If the supernovae rate was much 
higher in the past than at present, the plasma of the first formed objects must have been 
much more turbulent than that of presently observed, low redshift, star forming galactic 
molecular clouds. 

Turbulence generates stochastic magnetic fields (magnetic noise) at a faster rate than 
it does mean fields (Kulsrud & Anderson 1992). If the turbulence is strongly non-hehcal, 
the fields induced are confined to small scales (Kazantzcv 1968). However, if it is helical, 
induction of large scale magnetic correlations by the a effect occurs (Vainshtein & 
Kichatinov 1986). Astrophysical turbulence is mainly of a helical nature. Hence, we can 
expect that large scale correlations were induced by the high redshift, turbulent plasma. 

In this study, we explore the hypothesis that the magnetic fields observed in high 
redshift galaxies were created by small scale, stochastic, turbulent helical dynamos, rather 
than by mean field dynamos. We use a simple model of a gas cloud that is assumed to 
have collapsed at a high redshift 2;>10. At2;~10, the cloud would have had a low 
magnetization level and a high level of turbulence. Thus, it would have been similar to the 
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turbulent, low ionization, star forming molecular clouds observed in our galaxy, albeit with 
a much smaller initial magnetic field and a much higher turbulence level due to the higher 
star formation rate in the early universe. 

It is well known that shock waves produced by supernova explosions accelerate cosmic 
rays to energies > GeV (e.g., Blandford & Eichler 1987). Volk et al. (1989) showed that 
in all galaxies, the supernova rate is a direct measure of the cosmic ray intensity. We can, 
therefore, infer that cosmic rays were already present in considerable intensities in high 
redshift galaxies. We take into account phenomenologically, the efi:ect on turbulence of 
cosmic rays, supernova shocks, and powerful stellar winds from massive stars on turbulence 
by varying the turbulent parameters over a broad range, in order to take into account the 
uncertainies in our knowledge of high redshift structures. 

The linear evolution equations for the correlation function of magnetic fields for 
non-helical turbulence were derived nearly forty years ago by Kazantzev (1968). For helical 
turbulence, the corresponding equations were obtained twenty years later by Vainshtein 
& Kichatinov (1986). These equations are hnear in the magnetic correlations. Recently, 
Subramanian (1999) and Brandemburg & Subramanian (2000) derived the non-linear 
evolution equations for the magnetic correlations by taking into account the back-reaction 
of the Lorentz force on the plasma charges in the form of ambipolar drift. We solve the 
nonlinear helical evolution equations numerically for various values of the parameters that 
characterize the high redshift turbulent plasma. 

The paper is organized as follows. In section II, we give the evolution equations for the 
magnetic correlations. We describe the effects of the main parameters on the integration in 
section III . Finally, in section IV, we summarize and discuss our results. 
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2. MAGNETIC FIELD EVOLUTION EQUATIONS 



In this section, we summarize the nonhnear evolution equations for the magnetic 
field correlations (Subramanian 1999, Brandemburg & Subramanian 2000). The 
evolution equation for the magnetic field is given by the induction equation, d'B/dt ~ 
V X (v X B — 77V X B) , where B is the magnetic field, v the velocity of the fiuid, and 
rj is the Ohmic resistivity. The velocity v {— Vt + Vd) is the sum of an external 
stochastic field v^- and an ambipolar drift component v^, which describes the non-linear 
back-reaction of the Lorentz force. This back-reaction is due to the force that the 
ionized gas exerts on the neutral gas through collisions of the ions with the neutral 
atoms. It is assumed that vy is an isotropic, homogeneous, Gaussian random field with 
a zero mean value and a delta correlation function in time (Markovian approximation). 
Its two point correlation function is (v^ (x, i) (y, s)) = T^^{r)5{t — s), where 



T'^ (r) = Tjvjv (r) [5'^ - r'r^ /r'^ ]+Tll (r) (rV^r^) +C (r) e'^^rf (Monin & Yaglom 1975). 



The symbol () denotes ensemble averaging over the stochastic velocities, r = |x — y| , 

= x"^ — y\ Tll {f) and Tjvat {r) are the longitudinal and transverse correlation functions 
of the velocity field, respectively, and C (r) is the hehcal term of the velocity correlations. 
As the magnetic field grows, the Lorentz force acts on the fiuid. We assume that the 
fiuid responds instantaneously and develops an extra drift velocity, proportional to the 
instantaneous Lorentz force. We, thus, express the drift velocity as V£) = a [(V x B) x B] , 
where a — r/AiTpi, r the characteristic response time, and pi is the ion density. 

Consider a system whose size S ^ Lc, where Lc is the coherence scale of the 
turbulence, for which the mean field averaged over any scale is negligible. We take B to 
be a homogeneous, isotropic, Gaussian random field with a negligible mean average value. 
Thus, we take the equal time, two point correlation of the magnetic field as 




(1) 
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where 



3 



(2) 



(Subramanian 1999). The symbol () denotes a double ensemble average over both the 
stochastic velocity and B fields, Ml (r, t) and {r, t) are the longitudinal and transverse 
correlation functions, respectively, of the magnetic field, and H (r, t) is the helical term 

of the correlations. Graphically, can be represented as — > and Mjv, as 

t t . Hence, positive values of Ml and Mjv correspond to parallel vectors and negative 

values, to anti-parallel vectors. Since V • B = 0, we have Mm = (1/2 r) 9 (r^Mz,) / {dr) 
(Monin & Yaglom 1975). The induction equation can be converted into evolution equations 
for Ml and H : 

OMl 



dt 



dH 
'dt 



ir,t) 



or \ or 
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. d 
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- aN{r,t)ML{r,t)]], 



(3) 



(4) 



where 



and 



kn {r,t) =r] + TLL (0) - Tll (r) + 2aMi (0,t) , 
aN (r, t) = 2 C (0) - 2 C (r) - 4 a (0, t) , 



G (r) = -4 



r d 




\dr 


r 



+ ^7-[rTLL (r)] 
r^ dr 



(5) 
(6) 

(7) 



(Subramanian 1999). These equations form a closed set of nonlinear partial differential 
equations for the evolution of Ml and H, describing the evolution of magnetic correlations 
at small and large scales. The effective diffusion coefficient k,jv includes microscopic 
diffusion {r}), a scale-dependent turbulent diffusion [Tll (0) — Tll {'t)] , and ambipolar drift 
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2aMi (0, t) , which is proportional to the energy density of the fluctuating fields. Similarly, 
is a scale-dependent a effect, proportional to [2 C (0) — 2 C (r)]. The nonlinear decrement 
of the a effect due to ambipolar drift is iaH (0, t) , proportional to the mean helicity of the 
magnetic fluctuations. The G (r) term in equation (3) allows for rapid generation of small 
scale magnetic fluctuations due to velocity shear (Zel'dovich et al. 1983; Kazantzerv 1968). 
We are interested in the evolution of {r) since this function gives information about the 
coherence of the induced large scale magnetic field. A positive value of this function over 
a given length indicates that the field is coherent in this region. Therefore this length will 
be taken as the coherence scale of the induced field. Since is the correlation function 
of the tensor product of parallel vectors, evaluated at two points separated by a distance 
r, we can estimate the induced magnetic field intensity at all points where > as 
B^ML{r)/Ml/\0). 

3. TURBULENT STOCHASTIC DYNAMO ACTION 
IN HIGH REDSHIFT GALAXIES 

In order to study the evolution of the magnetic correlations due to the turbulent plasma 
in the high redshift objects, we integrated equations (3) and (4) numerically for different 
values of the parameters. We employed second order conservative finite differencing in 
space, with Neumann boundary conditions. In the time discretization, we used a second 
order Crank-Nicolson type method, except for the treatment of the non-linear terms. In 
these terms, we employed the values of Mi(0,t) and H{0,t) from the previous time-step, 
making the system of equations to be solved, linear in each time-step. The equations were 
solved by a few iterations of a relaxation procedure. The implicit treatment in the time 
discretization is important to avoid the severe stability constraints that would result from 
a fully explicit time discretization of the system. For the numerical results presented in 
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this paper, we employed a spatial grid with 5000 equally spaced grid-points. With this 
resolution, we were able to obtain convergence. Doubling the resolution led to graphically 
indistinguishable numerical results. 



3.1. Chciracterizing the High Redshift Plasmas 

We considered a cloud at ^ ~ 10 and followed the evolution of the magnetic correlations 
until z ~ 5 (~ 10^ years). The value taken for the cut-off scale of the turbulence, Zc ~ 1 AU, 
is similar to that for present objects (Zel'dovich et al. 1983). Assuming 3> /c, we 
studied the range of values 10 pc < Lc < 100 pc. We assumed that the height h of the 
turbulent eddies of the high redshift object is of the same order of magnitude as Lc- 
In order to estimate the correlation velocity Vc on the scale Lc, we used the expression 
Vj^ {yd Lc) ~ £, where e is the turbulent energy dissipated per unit mass per unit time. 
This expression assumes that the energy is dissipated on the order of a single rotation of 

1/3 

the eddies of size Lc at the angular frequency $1^14/ -^c- We then have Vc ~ {sLc) . 
Supernova explosions are a major contributor to the galactic turbulent energy. The energy 
associated with a supernova remnant in our galaxy is about 3 x 10^° erg, with about one 
third transformed into kinetic energy of the ambient gas. Larger values for the supernova 
remnant energy and the mass of the gas involved in the explosions, will produce higher 
turbulent velocities. We assumed that at rcdshifts 5-10, / explosions occurred every 5 
years and that the mass of the gas involved was 10 ^^M© (Zel'dovich et al. 1983). As noted 
above, the star formation and supernova rates were very high in the past. The indicated 
star formation rate from observations increased by a factor of ~ 50, in going from 2; ~ 
to 2; ~ 8 (see e.g., fig. 4 in Lanzetta et al. 2002). The expected values for / are then 
1 < / < 10. A value of / ~ 0.1 corresponds to the present supernova rate in our galaxy. 
We, thus, have £ ~ 0.3 x / cm^s"^. For the considered values of L^ the expected range 
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of values for Vc is 9.59 km s~^ i$ K i$ 96.5 km s"^. These values are 3-10 times larger 
than those in our galaxy (Zel'dovich et al. 1983). Assuming that the largest velocity 
corresponds to the largest eddy, we have fl ~ 10~^^s~^. We estimated that the baryon 
density is p„ (z) — p„ (0) (1 + ^ b, where p„ (0) is the present baryon density and 6 is a 
compression factor, which can be much greater than ~ 200 (virial collapse). In our galaxy, 
the particle density is ~ 1 cm~^ or p„ ~ 10~^^gcm~^. The average baryon density in the 
universe today is ~ 10~^°gcm~^. Thus, for our galaxy, the compression factor is 6 ~ 10^. 
We assumed that the cloud that we are studying in the interval 5 < 2; < 10, collapsed 
virially at a high redshift, creating a large b. Reasonable values for b are, then, in the range 
200 < 6 < 10 ^ Taking p„ (0) - 0.05 pc (0) , where pc (0) ~ 0.9 x lO'^^ g cm'^ is the present 
critical density (assuming a fiducial factor, h ~ 0.7, for the Hubble constant), we obtain 

4 X lO^^^gcm"^ < p^ (^z — 10) < 2.3 x 10~^^gcm~^ for the baryon density in our high 
redshift cloud. We estimated the ion mass density as pi ~ gpn, with 0.001 ^ g ^ I, which 
gives an ion density in the range 4 x 10~^^gcm~^ ^ Pi ^ 2.3 x 10~^^ gcm~^. 

At 2; ~ 10, the cosmic microwave radiation temperature was (1 + z) Tq ~ 30 K. For 

5 < 2 < 10, we considered plasma cloud temperatures in the interval 30K < T < lO^K. 
Using these values and estimating the thermal velocity of the ions as Vn — {?>kBT /nipY/'^ , 
we obtained lO^cms"^ ^ ^ lO^cms"-*^. Comparing these values with we see that we 
are dealing with mildly supersonic turbulence. 

Due to the the relatively low temperatures of the plasma, the ion-neutral collision 
cross section is CTj^ ~ 10~^^cm^ (NRL 2002). The ion- neutral collision frequency is 
^in — f^innnVth, giving 10~^^ < < 10~^°s~^. The electrical resistivity can be 
estimated as rj — (c^/^tt) (mg i^en/e^ rig) , where Ug is the electron number density, nie 
the electron mass, and Ugn = (o'enVe) Un is the electron-neutral collision frequency. 
Taking Ue — rii (charge neutrality), Tg ~ Tj, and using ~ {SkBTe/nie)^^'^ , we obtain 
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?7 ~ 5 X lO^cm^s""^, which is extremely small. The magnetic Reynolds number is 
Rm = LcVc/rj ~ 10^^ — 10^^, which means that at high redshifts, plasma turbulence 
was the main mechanism for diffusion and dissipation. Thus the first term in equation 
(5) can be neglected. Since the ion-neutral collision was the dominant interaction in the 
plasmas considered, we took the characteristic response time as r ~ u'^. The coefficient 
"a" in the non-linear terms in equations (5) and (6) can then assume values in the interval 
4.3 X 10 30 g-^cm^s < a < 2.5 x 10^^ g-^cm^s. 



3.2. Chciracterizing the Turbulence 



When studying low velocity (|vt| <^ velocity of sound) turbulence, it is usually 
assumed that the fluid is incompressible (V-Vt = 0). The functions T^n and Tll are, 
then, related in the way described by Subramanian (1999). When the above approximation 
is not valid (as is the case here), V x vt = is used and these functions are related by 
Tll — Tnn -\- f {dT^N I dr) (Monin & Yaglom 1975). The fluid flow correlation functions 
can be written as 



2C{r) 



h 



Tnn (r) = r > Lc, 



< r < Lc, 
lc<r < Lc, 



(8) 

(9) 
(10) 



with An = VcLc/3 (Vainshtein 1982). In our study, lc is much smaller than the numerical 
resolution used. We, therefore, considered Ml (0) = Ml {lc) ■ For free turbulence, we have 
p — 4/3, (the Richardson law) (McComb 1990). We take a range of values for p to take 
into account the uncertainties in the physics of high redshift objects 1 < p < 2. It is 
customary to take g = 2 for the helicity spectrum, but here, we shall allow for a more 
general dependence, using q > 1. 
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The required integration time can be estimated from the fact that when the kinetic 
energy density of the turbulence equals the magnetic energy density (i.e., Pn Vc /2 ~ Ml (Ic)), 
turbulence cannot supply more energy to create stronger magnetic fields. In the integrations 
that we performed, the growth saturated before this condition was reached. 

3.3. Discussion 

Both the size of the coherent region Lm and the induced intensity of the magnetic 
field Bm were studied. We estimated the value of Bm at all points where > as 

1 /2 

Bm {f) ~ Ml (r) /M^' (0) . We found that, in general, the magnetic correlations that result 
from the evolution of the turbulent kinematical dynamo in going from 2; = 10tO2; = 5, are 
independent of the initial field correlations. 

We investigated the following sets of turbulent parameters: 1) Lc — 33 pc and 
Vc — 45kms~^; and 2) Lc = 81 pc and Vc — 96kms~^. In both cases, we took p — 4/3, 
g = 2, ~ 5 X 10'^ cm^ s~^, and a ~ 9.76 x 10^^ cm^ sg^^. For the first set of parameters, we 
obtained Bm ~ 1-1 x 10~^ G and Lm — 1.7 kpc and for the second case, Bm ~ 1-4 x 10~^ G 
and Lm ~ 5.4 kpc. 

In general, Bm at 2; = 5 is sensitive to the values of a and Vc- The parameter a depends 
on both Pi and such that it is not possible to discriminate the dependence of our results 
on each of these factors, independently. The length Lm depends mainly on the value of V^. 
Different values of p and q change Bm and Lm only slightly. 

In Figures (l)-(4), we show (r) at 2; = 5 for different values of the parameters. 

In Figure 1, we plotted Ml (r) as a function of r for p = 1.1, 1.33, and 1.66. Using 
Lc = 81 pc, Vc = 96kms~^, a = 9.76 x lO^^cm^sg"^, and g = 2, we see that Bm and Lm 
are somewhat larger for smaller values of p. 



2000 4000 6000 8000 10000 

r (pc) 



Fig. 1. — Final value of M/^(G^) as a function of r(pc) for p = 1.11, 1.33, and 1.66; q = 2; 
Lc — 81 pc; Vc — 98kms~^; and a — 9.76 x 10^® cm^ sg~^. We see that Bm and Lm are 
somewhat larger for smaller values of p. 
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In Figure 2, we show Ml (r) as a function of r for g = 1.8, 2, and 2.5, using p — 1.33 
and the same values for Lc, Vc, and a as in Figure 1. The values of Bm and Lm are a little 
larger for smaller values of q. 

In Figure 3, we plotted Ml (r) as a function of r for a = 3.2 x 10^^ cm^sg"^, 9.7 x 
10"^^ cm^sg"^, and 2.43 x 10"^° cm^sg"^; p = 1.33; and q = 2. We used the same values for 
Lc and Vc as in Fig. 1. We see that the smaller the value of a (high ion density and/or 
high ion-neutral collision frequency), the larger the value of Bm- However, Lm is almost 
insensitive to the value of a. 

In Figure 4, we plotted Ml (r) as a function of r for g = 2, p = 1.33, Lc — 81 pc, 
a = 3.2 X lO^^cm^sg"^, and Vc = 30kms~^, 45kms~^, and 96kms~^. We see that large 
values of these two parameters produce high values of Ml as well as large coherence lengths. 

Finally, in Figure 5, we plotted the evolution of Ml as a function of t at ro — 112 pc, 
going from z = 10 to z = 5, for Ml {vq, i = 0) = 1.5 x 10-=^^ G^, 1.5 x 10"^^ G^, and 
1.5 X 10~^^ G^. We see that after a short period of time, f ~ 5 x 10^ years. Ml reaches its 
saturation value, M^^^ ~ 10~^ G, independent of its initial value. 

Throughout the integration time, the kinetic energy density of the fluid was greater 
than that of the magnetic energy. The magnetic energy density is given by Eb = B'^/Stt, 
which for B ~ 10~^ G (see Fig. 5), has a value of Eb ^ 10~^^ ergcm"^. The kinetic 
energy density is given by Ey — p^^c^ 2. For Vc — 98 kms"^ and the range of values used 

for p„, the kinetic energy, 3 x 10~^^ ergcm"^ < Ey < 10^^ ergcm^^, was many orders of 
magnitude greater than the magnetic energy. 
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Fig. 2. — Final Mj^(G^) as a function of r(pc) for q = 1.5,2, and 2.5; p = 1.333; Lc = 81 pc; 
Vc — 98kms; a — 9.76 x lO^^cm^sg"^; and p — 1.33. For smaller values of g, Bm and Lm 
are somewhat larger. 



-16- 



1.2e-10 
le-10 



9. 6e-ll - 



4e-ll 







1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
— 


1 1 1 1 1 1 1 1 1 


1 1 


\ 


a = 3.2e+39 

a = 9.7e+39 

a = 2.43e+40 




N 






N 
N 

N 

S 

— N 






N 

\ 

..^ ^ 

_ ^ ■ ^ . \ 










-1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




—1— J — 



1000 2000 3000 4000 

r(pc) 



5000 6000 



7000 



Fig. 3. — Final values of Ml{G^) as a function of r(pc) for g = 2, p = 1.333, Lc — 81 pc, and 
Vc = 98kms~^ and three different values for a : 3.2 x lO^^cm^sg"^ 9.76 x lO^^cm^sg"^ 
and 2.43 x 10^^ cm^ sg~^. For smaller values of a, Bm is larger, while Lm remains practically 
unchanged. 
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Fig. 4. — Final values of Mi(G^) as a function of r(pc) for q = 2, p = 1.33, L,. = 81 pc, 
and a — 9.76 x lO^^cm^sg"^ and three different values for V^:30kms~^, 45kms~^, and 
96kms~^. We see that for larger values of Bm and Lm are larger. 
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Fig. 5. — Values of MiiC^) as a function of t(years) at ro = 112 pc for three different initial 
values: Ml (ro, i = 0) = 1.5 x IQ-^^G^, 1.5 x 10"^^ G^ and 1.5 x IQ-^^ q2 ^g^^ ^ ^ ^ -^^^^ 
q = 2, Lc = 81 pc, Vc = 98kms, and a = 9.76 x 10^^ cm^ sg~^. Independent of the initial 
values of M^, saturation occurs at a time i ~ 5 x 10^ years, when ~ 10~^ G. 
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4. CONCLUSIONS 

In this work, we studied the problem of the origin of strong coherent large-scale 
magnetic fields, observed in low-redshift galaxies and previously thought to have been 
created by the mean field dynamo. Since doubts have been cast on the mean field dynamo 
as the source of these fields due to their observation in high redshift objects, where the 
dynamo would not have had suficient time to operate, we investigated here the stochastic 
helical dynamo as such a source. We showed that these fields can be generated in the 
plasmas found in very high redshift objects. 

The generation of strong large scale coherent fields is a highly non-linear magnetohy- 
drodynamical problem, which depends upon many factors. Here we discussed a possible 
mechanism for the generation of large scale fields, namely the non-linear evolution and 
diffusion of magnetic noise. Magnetic noise becomes coherent on a scale which is larger than 
that of the turbulence due to the presence of non linear terms in the evolution equations for 
the magnetic correlations. 

We found that for realistic turbulent parameters, it is possible to generate the magnetic 
fields, observed at high redshifts . Between 2; ~ 10 and 2; ~ 5, the primordial plasma 
is strongly turbulent and partially ionized The magnetic field intensities reached are 
independent of initial correlations. 

We considered a very simple model for the generation and evolution of the magnetic 
correlations. The dependence of the resulting magnetic field intensity on the charge 
composition of the primordial plasma, suggest that the reionization and star formation 
processes played an important role in determining the features of the magnetic fields 
detected in high redshift objects. In a forthcoming work we shall address the evolution 
of the magnetic correlations, considering a time dependent ion density as well as other 
nonlinear processes, such as the the Hall effect. Other turbulent scenarios, in addition to 



-20- 



the homogeneous and isotropic one considered here, will be treated as well. 
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